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Extending previous work [T. Pusztai, L. Ratkai, A. Szallas, and L. Granasy, Phys. Rev. E 87 , 032402 
(2013)], we have studied the formation of eutectic dendrites in a model ternary system within the frame¬ 
work of the phase-field theory. We have mapped out the domain in which two-phase dendritic structures 
grow. With increasing pulling velocity, the following sequence of growth morphologies is observed: flat front 
lamellae —> eutectic colonies —> eutectic dendrites —> dendrites with target pattern —> partitionless dendrites 
—> partitionless flat front. We confirm that the two-phase and one-phase dendrites have similar forms, and 
display a similar scaling of the dendrite tip radius with the interface free energy. It is also found that the 
possible eutectic patterns include the target pattern, and single- and multiarm spirals, of which the thermal 
fluctuations choose. The most probable number of spiral arms increases with increasing tip radius and with 
decreasing kinetic anisotropy. Our numerical simulations confirm that in agreement with the assumptions of a 
recent analysis of two-phase dendrites [S. Akamatsu, S. Bottin-Rousseau, G. Faivre, and E. A. Brener, Phys. 
Rev. Lett. 112, 105502 (2014)], the Jackson-Hunt scaling of the eutectic wavelength with pulling velocity is 
obeyed in the parameter domain explored, and that the natural eutectic wavelength is proportional to the 
tip radius of the two-phase dendrites. Finally, we find that it is very difficult/virtually impossible to form 
spiraling two-phase dendrites without anisotropy, an observation that seems to contradict the expectations 
of Akamatsu et al.. Yet, it cannot be excluded, that in isotropic systems two-phase dendrites are rare events 
difficult to observe in simulations. 

PACS numbers: 68.70.+W, 81.10.Aj, 81.30.Fb 

I. INTRODUCTION 

A broad variety of systems show multi-arm spiral¬ 
ing, including spiral galaxief^ banded spherulites in 
polymeric systemiP, biological excitable system^ 3 ’® os¬ 
cillating chemical reactiond^l, spiraling growth edges 
on flat crystalline surfaced, ridges in sputtered high- 
temperature superconducting thin films 8 certain semi¬ 
conductor mate rials g rown by molecular beam epitaxy 10 , 
binary eutectics 1 ®^, an d helical Liesegang systems®^, 
and the recently discovered ternary eutectic dendrites 
(Fig. [l]fSI. Although the formation mechanism of spi¬ 
ral structures excite the fantasy of scientists for a long 
time, a general explanation is not available partly due 
to the diversity of the underlying physical phenomena. 

While the details differ in various realizations of spiral 
growth, diffusion and phase separation often play a role 
in the respective models. For example, the aggregation 
of starving cells is controlled by propagating spiral waves 
of a chemo-attractant, often yielding multiarmed spiral 
patterns in the case of slime mold 15 . In binary eutectics, 
spiraling has been associated with a specific anisotropy of 
the free energy of the solid-solid interfaced, the presence 
of screw dislocation^®®, or recently to fingering driven by 
osmotic fWd In contrast, the newly discovered spi¬ 
raling ternary eutectic dendrites emerge from the inter- 
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play of two-phase solidification with the Mullins-Sekerka- 
type diffusional instability caused by a third component, 
whose solubility differs in the solid and liquid phases®®. 
The complex microstructures observed in some of the 
ternary alloys®^ are also suspected to emerge due to the 
presence of two-phase dendrited. It is worth noting that 
besides their scientific interest and pleasing view the lat¬ 
ter class of spiraling / helical structures have been identi¬ 
fied as possible means of creating chiral metamaterials for 
optical applications via eutectic self-organization 20 " 23 . 

Models of spiral growth range from wave theorj®EH, 



FIG. 1. Spiraling eutectic dendrites during bulk-sample 
directional solidification of a SCN-DC-NA ternary eutectic 
alloy. Growth rates: (a) 0.28 /ims -1 ; (b) 0.95 /ims -1 . The 
bar corresponds to 20 pm. (Reprinted with permission from 
RefP^© 2010 American Physical Society.) 
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via the FitzHugh-Nagumo (FN) theor y for excitable 
medicP' 5 and reaction-diffusion model ® * 26 * 27 ! to the 
Ginzburg-Landau/phase-field type model®®. Studies 
of biological excitable media, relying on the FN model, 
have clarified essential features of spiral growth® The 
mechanism by which the multi-arm spirals do form is 
the attraction of single spirals rotating in the same di¬ 
rection, whereas the number of spiral arms is associ¬ 
ated with the ratio of the single spiral period to the 
refractoriness of the medium. Apparently, in such sys¬ 
tems the multi-arm spiral structures are unstable un¬ 
less the excitability is sufficiently low. Some other mod¬ 
els provide single spiral structures exclusively * 28 * 29 *. Re¬ 
markably, even pure confined systems were found to dis¬ 
play chiral symmetry breaking, and thus spiral growth 
in phase-field simulation®! A recent work based on 
Cahn-Hilliard type reaction-diffusion model of helical 
Liesegang systems indicate that single and multiple he¬ 
lixes may occur and the fluctuations choose from the pos¬ 
sible configurations!®. A similar scenario has been out¬ 
lined in a recent work of us in which ternary eutectic 
dendrites may realize a variety of eutectic patterns in¬ 
cluding a target pattern and multiple spiraling motifs 
of which thermal fluctuations of the system choos®! 
In our work, we have been looking for a suitable phys¬ 
ical model and shown that a simple ternary extension 
of the standard binary phase-field theory is able to de¬ 
scribe the formation of eutectic dendrites, offering thus 
the first model that is able to capture details of such ex¬ 
otic growth patterns. We note, furthermore, that this 
is another remarkable success of the phase-field theory: 
Without any specific adjustments to make spiraling eu¬ 
tectic dendrites, a minimal and fairly standard version 
of the phase-field theory of ternary solidification proved 
capable of recovering these rather complex solidification 
morphologies. We have shown that (a) in many respects 
eutectic dendrites behave like the single-phase dendrites 
(e.g., shape, tip radius vs. interface free energy, etc.); (b) 
there is a multiplicity of eutectic patterns, which lead to 
steady-state dendritic growth under nominally the same 
conditions; these are the target pattern, and single- to 
multi-arm spiraling structures; (c) it is the fluctuations 
that choose from these competing growth modes; (d) the 
number of spiral arms tends to increase with increasing 
dendrite tip radius, however, with some stochastic scat¬ 
tering. In a more recent analytic study, Akamatsu and 
coworkers 31 have addressed dendrite formation in ternary 
eutectic systems. They argued that assuming (i) A ~ R 
and (ii) X 2 v = const, (the latter is the Jackson-Hunt re¬ 
lationship of eutectic solidificatioi®®, where A, R 1 and 
v are the eutectic wavelength, the tip radius of the den¬ 
drite, and the tip velocity, respectively), the dendrite se¬ 
lection mechanism differs from the one seen in binary 
systems suggesting that in the case of ternary spiraling 
dendrites anisotropy is not neccessary for steady-state 
dendritic growth. Remarkably, assumption (ii) seems to 
be in contradiction with the analytic predictions by Liu 
et aZ.® for higher concentrations (5%) of the third com¬ 


ponent, obtained assuming a rotational paraboloid shape 
for the two-phase dendrite. However, the approach of Liu 
et al. appears to contain questionable assumptions that 
the lamella spacing is much smaller than the tip radius 
(expected true rather in the case of eutectic colonie®!), 
and the neglection of the effect of curvature of the den¬ 
drite surface on eutectic solidification. It is yet desirable 
to clarify these contradictory issues. 

Herein, we focus our attention to the ternary eutectic 
dendrites and extend our previous phase-field study 30 in 
the following directions: (a) a detailed investigation of 
the stochastic eutectic pattern selection mechanism; (b) 
testing of scaling laws (i) and (ii); and (c) interaction of 
neighboring eutectic dendrites. The structure of our pa¬ 
per is as follows: In Section II, we give a brief description 
of our model and justify the model parameters used dur¬ 
ing the simulations. In Section III, we review details of 
the performed numerical simulations, whereas in Section 
IV we present our results for the formation of two-phase 
dendrites, the scaling laws (i) and (ii), the pattern se¬ 
lection mechanism, the role played by the anisotropies, 
and the dendrite-dendrite interaction. Finally, we sum¬ 
marize our results in Section V and offer a few concluding 
remarks. 


II. MODEL DESCRIPTION 

We adopt here a simple ternary phase-field model 
used in our previous study®!. It ca n be obtained 
as a straightforward generalization 30 ® of the phase- 
field model of binary solidificatioii 37 !® whose variants 
and extensions h ave been used successfully? for describ¬ 
ing homogeneous!^ 1 ® and heterogeneou®*® nucleation, 
and cr ystal gr owth including dendritic® and eutectic 
freezing®®®, and even polycrystalline solidi ficat ion 
treating nucleation and growth on equal footing®®®. 
Such models turned ou t to be quantitative on the 
nanoscale for nucleatioi l 41 * 42 * 45 *, whereas qua litative re¬ 
sults were obtained for broader interface®®®. A sim¬ 
ilar model has been used to describe the formation of 
eutectic colonies in ternary system®! 

Herein, the ternary model is solved in a dimension¬ 
less form; accordingly, the size scale is determined by the 
interface thickness chosen for redimensionalizing the re¬ 
sults. Considering that the physical interface thickness 
is about 1-2 nm, the present results refer to nano-scale 
eutectic structures, expected to form at extreme under¬ 
coolings, which can be realized by e.g., via crystallizing 
amorphous alloys (see e.g.,®®). In contrast, the exper¬ 
iments refer to small undercoolings, where other models 
may be considered more appropriate (e.g.,®®). 

The free energy functional has been obtained as a 
simple generalization of the standard binary phase-field 
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FIG. 2. The free energy functions, f a ,i( c), at a tempera¬ 
ture displayed above the Gibbs simplex: (a) bulk solid; (b) 
bulk liquid; (c) liquid and solid free energy functions with the 
tangential plane which specifies the equilibrium solid and liq¬ 
uid phases; (d) Yellow circles stand for the three equilibrium 
phases, the ci- and C 2 -rich solid phases and the liquid phase. 
The white circle indicates the composition of the initial liquid 
typically used in the simulations. 


model (see e.g. RefsP^D): 

F[i>,c] =/ ~2 ( V "F w 9i ( f > ) + (1 -J#))/;( c )+ 

+P(4>) (/ s (c) + e f £i=i (Vg) 2 ) dV. (1) 


Here </>(r, t ) £ [0,1] is the phase held, which is 0 for bulk 
liquid and 1 for solid, c = ( 01 , 02 , 03 ) are the concen¬ 
tration fields, the quartic function g(<p) = cj> 2 ( 1 — 0) 2 /4 
ensures the double-well form of F, while function p((j>) = 
0 3 (10 — 15 4> + 6 4> 2 ) switches on and off the free energy 
densities of the solid and liquid phases. We choose ideal 
solution model for the bulk liquid phase [Fig. [2](b)]: 

3 

fi (c) = Ci [ fi + lo § c i] > ( 2 ) 

i =1 


where f\ = 0 , and the regular solution model for the solid 
[Fig- da)]: 

3 1 

fs(c) =5> [fi + log Ci] + - SlijCiCj. (3) 

i=1 


The equations of motion (EOMs) derived variationally 
have the following forms for the phase held 


= e(n )M (j> - wg' (<j>)+ 

+p'mfi(c) ~ / s (c))-p'W| £-=i(Vg) 2 ,1 (4) 


and for the concentration helds 





( 5 ) 


where SF/Sci is the functional derivative of the free en¬ 
ergy with respect to concentration held Cj. Here the 
£ ; Ci = 1 constraint is automatically satished by our 
choice of 1 and —0.5 for the diagonal and off-diagonal 
elements of the 3x3 mobility matrix, M c . Note that 
diffusion is switched off in the bulk solid. Owing to this 
and the ideal solution thermodynamics assumed in the 
bulk liquid, phase separation is possible in only the solid- 
liquid interface layer. 

In a few test cases that explore the effect of concentra¬ 
tion fluctuations, we included a conservative (hux) noise 
for the concentration helds, via adding VEj to the RHS of 
the respective EOMs [Eqs. ([5])], where the vector Ej(r, t) 
is a random current of component i. Since local mass 
conservation (£b d = 1 ) has to be obeyed even in the 
presence of fluctuations, the sum of the divergencies of 
these random currents must be zero, i.e., £, VEj = 0 
has to be satished. Local mass conservation is realized 
by the construction E; = where represents 

that part of the total current Ej, where the flow of com¬ 
ponent i is compensated by the back-flow of component 
j as dictated by the linked term £y(r, t) = —£ji(r,t) in 
the respective expression for E ? . When solving the dis¬ 
cretized equations, the random currents ^ = —£ji were 
realized by exchanging random amounts of components i 
and j between each possible pairs of neighboring cells and 
for all the i-j combinations. Random numbers of Gaus¬ 
sian distribution and zero mean were used to determine 
the amount of the components exchanged. 

We have considered in this work a kinetic anisotropy of 
cubic symmetry, represented by the following orientation 
dependent coefficient in Eq. (4): 


e(n) 


(1 - 3e 4 ) 



4e 4 

1 — 3e 4 



( 6 ) 


Here n = V^/|V<(>|, whereas parameter e 4 determines the 
magnitude of anisotropy. As this anisotropy acts at the 
solid-liquid interface independently of composition, the 
same kinetic anisotropy applies for the a-liquid and /3- 
liquid interfaces, as well as to the two-phase-solid liquid 
interface. Owing to neglecting the anisotropy of the solid- 
liquid interface free energy, the present computations are 
regarded as qualitative. 


III. NUMERICAL SIMULATIONS 

The equations of motion were solved numerically in 
a dimensionless form using the finite difference method 
and explicit time stepping. The simulations were per¬ 
formed in a massively parallel manner partly on a com¬ 
puter cluster of 928 CPU cores and on several graphics 
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processing unit (GPU) cards housed at the Wigner Re¬ 
search Centre for Physics, Budapest, Hungary. In a few 
cases, where extreme large simulations were required, we 
used the supercomputer Hitachi Super Technical Server 
SR16000 Model Ml of the Center for Computational Ma¬ 
terials Science, at the Institute of Materials Research, 
Tohoku University, Sendai, Japan. 

The simulations have been carried out under condi¬ 
tions that correspond to directional solidification. A 
temperature gradient was imposed via making the free 
energy of the solid dependent on position: f s ,i,z = 
/ s ° — z (df St i/dz). Here z is the spatial coordinate along 
the direction of the pulling velocity v p . In order to model 
the pulling of the sample, the contents of the arrays (f> and 
Cj have been shifted by one voxel back in the z direction 
in each [(dx / dt) / v p ] th time step, while prescribing the 
boundary conditions <fi = 0 and c = Co on the high T 
and no-flux boundary conditions on the low T side of the 
simulation box. 

To make simulations that are sufficiently large in the 
direction of pulling (5) more efficient, only Eqs. ([5]) have 
been solved far ahead of the solidification front, where 
the phase field is small (< 10 -8 ) to be regarded as bulk 
liquid. This approach allowed us to handle essentially 
the whole diffusion zone ahead of the solidification front, 
which then ensured that the dendrite morphology became 
independent of the length of the simulation box (N z ). 

As in our previous studj®, the majority of our sim¬ 
ulations has been started by including a slab of solid of 
length (2/9 )L Z into the simulation box (L z = N z dx is the 
length of the simulation box in direction z) with a small 
hump at the center, whereas the initial composition of the 
solid has been set to (ci) = ( 02 ) = 0.455 and C 3 = 0.09 
[see Fig. id)], realized by a random transversal (x — y 
plane) distribution of the two solid phases of typical size 
scale close to the natural wavelength of eutectic growth. 
The random transversal distribution represents here the 
effect of all the thermal fluctuations from the nucleation 
of the eutectic grain to the starting of the simulation. We 
have opted for this initial condition, because simulations 
that follow the formation of the two-phase dendrites from 
fluctuation-induced emergence of nuclei are prohibitively 
time consuming. Evidently, this might be a crude ap¬ 
proximation. Therefore, we tested in a few cases how 
close one can get to the results obtained from this manner 
if we start the solidification from a small but supercriti¬ 
cal (larger than the critical size for nucleation) one-phase 



FIG. 4. Compositions at which two-phase dendrites of target 
or spiral patterns have been observed at v p = 0.2 (red dots). 


solid sphere, while adding flux-noise to the EOMs of the 
concentration fields as described in RefP^. As will be 
seen, the probabilities of formation of the individual eu¬ 
tectic patterns obtained via these two routes are fairly 
close, indicating physical consistency. 

Reference conditions: These were used unless stated 
otherwise. Time and spatial steps: dt = 0.0025 and 
dx = 1.0. Grid: N x x N y x N z = 96 x 96 x 612. Size of 
the simulation box: Li — Ni ■ dx , where * = x,y, or z. 
Composition: Ci = C 2 = 0.455, and C 3 = 0.09. Parame¬ 
ters of free energy densities: fi ;i = 0; f^°- = —0.9640 
and df s j/dz = 3.677 x 10 ” 4 . At t = 9000 (steady- 
state), this yields = —0.9140 at the dendrite tip, 
corresponding to a rather substantial relative undercool¬ 
ing of A T r = (Tl — T)/(Tl — Tg) S 3 0.78, where Tg and 
Tg are the dimensionless liquidus and solidus temper¬ 
atures corresponding to the conditions at the dendrite 
tip. Pulling velocity: v p = 0.2, whereas f2i >2 = 3.05, 
^ 2,3 = ^ 3,1 = 3.0. M 0 = 1.0. Solid-liquid interface free 
energy (isotropic): 7 sx,o = 0.0147. Kinetic anisotropy: 
e 4 = 0.3. Coefficients of the square-gradient terms: 

= 0.75, ej? = 0.4, and the free energy scale, w = 0.0469. 
Unless stated otherwise no noise has been added to the 
EOMs. 



FIG. 3. Lateral view of the simulation box. Grid size: N x x 
N y x N z = 96 x 96 x 612, unless stated otherwise. The dashed 
line indicates the temperature distribution. The yellow and 
red domains stand for the solid phases, while the liquid is 
made transparent. 


IV. RESULTS AND DISCUSSION 

In a previous work on spiraling eutectic dendrite^^, 
we explored the parameter space defined by composition, 
temperature gradient, pulling velocity, interfacial free en¬ 
ergy, and kinetic anisotropy, and optimized the condi¬ 
tions for growing two-phase steady-state dendritic struc¬ 
tures. Herein, we extend these studies to significantly 
broader parameter ranges. We examine, furthermore, 
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FIG. 5. Solidification morphology and pattern formation as a function of dimensionless pulling velocity v p . The ci and C 2 -rich 
solid solutions are colored red and yellow, respectively, whereas the liquid is transparent, and purple stands for ci « C 2 - The 
front view (top row), the longitudinal section (central row) and transverse cross sections (bottom row) are displayed. With 
increasing velocity, the sequence of growth morphologies is: flat front lamellae —» eutectic colonies —» eutectic dendrites —> 
dendrites with target pattern —> partitionless dendrites —¥ partitionless flat front. At extremely low and high velocities planar 
fronts develop. Above v p fts 0.35 solidification takes place without partitioning. At extreme high velocities ( v p > 0.8) full solute 
trapping occurs. Remarkable is the gradual transition from the (usual) lamellae parallel to the temperature gradient seen at 
low velocities to lamellae perpendicular to the temperature gradient seen at v p 0.3. 


the relationship between the tip radius and the eutec¬ 
tic wavelength, inve stigat e the validity of Jackson-Hunt 
scaling ( X 2 v = const IMM), present a detailed character¬ 
ization of the morphology of the forming eutectic struc¬ 
tures, including the 3D geometry of the dendrite tip, and 
the internal distribution of the constituent phases, and 
investigate whether anisotropy is indeed not needed for 
the appearance of two-phase dendrites. 



FIG. 6. The wavelength of eutectic solidification as a 
function of pulling velocity. (Symbols of different colors 
denote series of simulations with different initial patterns.) 
Note the reasonable agreement with the Jackson-Hunt scal¬ 
ing, A oc v p 1,7 2 (solid line). 


A. Domain of ordered two-phase dendrites 

To identify the region, in which spiraling eutectic den¬ 
drites do form in the phase diagram, we have investi¬ 
gated the solidification morphology in the neighborhood 
of the operating point [c = (0.455, 0.455, 0.09) and 
v p = 0.2], reported in Ref. 30 . The respective domain 
is shown in Fig. [4] which indicates that spiraling eutectic 
dendrites may be observed in an oval region that extends 
to 0.35 < Ci < 0.56 along the line C3 = 0.09, whereas 
it covers the range 0.05 < C3 < 0.17 along the line of 
symmetric compositions (ci = 02). 

B. Jackson-Hunt scaling 

To see whether the \ 2 v = const, relationship applies 
indeed, we have investigated the solidification morphol¬ 
ogy as a function of pulling velocity at the composi¬ 
tion c = (0.455, 0.455, 0.09). For low pulling velocities 
(v p < 0.01) one finds a planar front with a disordered 
lamellar pattern (see Fig. |5j|. Increasing v p resulted in a 
dendritic morphology yet with disordered lamellar pat¬ 
tern. Around v p ~ 0.2 ordered eutectic patterns form 
such as the target and spiraling patterns. At higher 
v p ~ 0.3, solidification takes place via forming lamellae 
perpendicular to the temperature gradient. Above this 
velocity partitionless dendrite forms indicating that so¬ 
lute trapping does not occur for C3 at the same level as for 
the other two components. At even higher pulling veloc- 
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FIG. 7. Eutectic wavelength vs. tip radius, (a) Velocity 
is varied, while other parameters are fixed. Different sym¬ 
bols stand for runs with different initial patterns (different 
initializations of the random number generator), (b) Varying 
the interface free energy 'ysLfi or the kinetic anisotropy £4 
while keeping the velocity constant. (Here circles, triangles 
and pentagons indicate the number of spiral arms observed in 
simulations where 'ysLfi was varied, while diamonds stand for 
results from simulations where the kinetic anisotropy 64 was 
varied. In the experiments ysn.o and £4 are fixed, under such 
conditions the two characteristic lengths are roughly propor¬ 
tional to each other, A oc Ru p , as assumed in Ref.™ However, 
the tip radius can be influenced by varying either the mag¬ 
nitude of the interfacial free energy or the kinetic anisotropy. 
The dashed lines are to guide the eye. 


ities ( v p ~ 2.0) partitionless growth with a flat interface 
has been observed. The results shown in Fig. [5] suggest 
that the transition from flat lamellar eutectic structure 
to the partitionless growth with flat interface with in¬ 
creasing velocity happens via the following sequence: flat 
lamellar front —> eutectic colonies —> eutectic dendrites —> 
dendrites with target pattern —» partitionless dendrites 
—I partitionless flat front. Remarkably, during the tran¬ 
sition from a flat surface with lamellae perpendicular to 
it (seen at small velocities) to a flat surface at full solute 
trapping (at high velocities), happens via a stage preced¬ 
ing solute trapping, in which a nearly ID eutectic pattern 
evolves: lamellae that are perpendicular to the tempera¬ 


ture gradient. Similar transition (lamellae parallel with 
the temperature gradient to lamellae perpendicular to it) 
has been observed w hen inc reasing the velocity in binary 
eutectic simulation d 38|40 | 58 J . 

Apparently, the natural length scale of eutectic solid¬ 
ification follows the Jackson-Hunt scaling (c/. symbols 
and solid line in Fig. [6]), confirming thus the assumption 
made in Ref^S. The other assumption of RefPS that the 
tip radius and the eutectic wavelength are proportional 
A oc R t ip is satisfied likewise provided that the other pa¬ 
rameters are kept constant [see Fig. [7](a)] . In our case, 
the proportionality coefficient A = (0.26 ± 4.5 %)R t i P is 
much smaller than in the experiments A ss 1.337?4, r pH. 
One can however, tune the the ratio via changing the 
interface free energy or the kinetic anisotropy: If the ve¬ 
locity is kept constant and either the interface free en¬ 
ergy is increased or the kinetic anisotropy decreased, the 
tip radius increases, whereas the eutectic wavelength re¬ 
mains roughly constant [Fig. 0b)P. Accordingly, one 
could move towards the experimental conditions via de¬ 
creasing the interface free energy or increasing the kinetic 






FIG. 8 . Two-phase two-armed spiraling dendrite grown un¬ 
der the reference conditions defined in Section III: (a) spi¬ 
raling motif on the surface; (b) the helical structure formed 
by one of the solid phases, (c) contour lines showing the 
transverse sections at lOdcc distances; (d) longitudinal sec¬ 
tion (dots), the best fit parabola (dashed), and the curve 
z = Zmax — \x\ v fitted to it (solid line); (e) maximum of C3 at 
the tip vs time. 
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anisotropy. Presumably, the presence of anisotropic in¬ 
terface free energy could also be used to tune X/R t i P , 
however, its investigation is left for the future. 


C. Dendrite morphology 


A typical two-phase dendrite is shown in Figs. Ufa) 
and[8]jb). It has a rounded square-like transverse section 
in the x — y plane [Fig. [8jc)] ; whereas in the fin direc¬ 
tions (e.g., x — z), the longitudinal profile can be fitted 
reasonably well by z = z max — \x\ v , where z max is the 
tip position, x is the distance from the axis of the den¬ 
drite, and v is ~ 1.58 ± 0.05 [Fig. [8j[d)] . This exponent 
is somewhat lower than v = 5/3 predicted theoretically 
for single-phase dendrites®®, or the v = 1.67 obtained ex¬ 
perimentally for xenon dendrites®®®!]. (j n our wor ] S) the 
perimeter of the dendrite has been defined as the contour 
line <f> = 0.5.) The quasi steady growth form (not a true 
steady-state as the eutectic pattern rotates in the coor¬ 
dinate system moving with the dendrite tip) has been 
established after a transient characterized by decaying 
oscillations of the tip-radius, tip-temperature, and the 
maximum of C3 at the tip, as shown in Fig. [8](e). These 
findings are in a reasonable agreement with a detailed 
characterization of tip shapes of single-phase dendrites 
by phase-field simulationJ^il, which indicates that there 
is a crossover between a corrected paraboloid valid close 
to the tip and the xenon-like behavior on the tail; there¬ 
fore, the results depend (here too) somewhat on the range 
of distances one choses to fit the shape. 

To see how far the two-phase dendrites follow the be¬ 
havior of single-phase dendrites, we have varied the mag¬ 
nitude of the solid-liquid interface free energy (jsl) via 
changing the free energy of the single component solid- 
liquid interface ( ; ysL,o )• After the dendrite reached its 
quasi steady-state, we evaluated and averaged the tip ra¬ 
dius in the fin directions ( Ru P )• The results are summa¬ 
rized in Fig. [9I It was found that R tip oc 7s|° 0 ±H ' 01 , which 

is very close to the theoretical relationship R t i P oc 
derived for single-phase dendrites (see e.g.p*® 01®. It 
probably indicates that the chemical contribution is neg¬ 
ligible to the interfacial free energy 7 sl- We should also 
call attention to the fact that the shape of the dendrite 
(and the tip radius) is independent of the two-phase pat¬ 
tern forming on the surface of the dendrite: The target 
pattern, and single- and multiple spiraling motifs appear 
to coexist under the same conditions (for further details 
see Sub-section IV.D below), yet the probability of find¬ 
ing a larger number of V arm increases with increasing tip 
radius. 

Another way to influence the shape of the dendrites 
is to change the kinetic anisotropy (£4). The results 
are shown in Fig. 10 a): decreasing anisotropy yields an 


increasing tip radius, a change apparently initiating a 
larger number of the spiral arms. The variation of the ki¬ 
netic anisotropy influences the shape of the cross-section 
in the direction of the fins. The exponent v describing 


100 




10 


— fit 

o 0 

X 1 

0 
A 

4 

O 5 






m = 0.50 ±0.01 
b - 4.16 ±0.02 


0.01 


0.1 


7 sl,o 


FIG. 9. Tip radius vs. solid-liquid interface energy. Results 
obtained with the same initialization but at different Tsn.o) 
values are depicted by different empty symbols indicating the 
actual eutectic pattern. In contrast, the black symbols stand 
for 20 additional simulations, which were performed under the 
same condition as the underlying empty ones, except that the 
initial random spatial distribution of the solid phases was dif¬ 
ferent (achieved by usig different initializations for the random 
number generator). These simulations realize various eutectic 
patterns, yet they fall on the same master curve, indicating 
that the shape of the dendrite is essentially independent of 
the eutectic pattern. 


the shape of the dendrite tip changes from ~ 1.5 ±0.1 to 
~ 2.1 ± 0.1 with decreasing anisotropy [see Fig. [U: b )]; 
i.e., it varies roughly between the experimental value for 
xenon (1.67) and the rotational paraboloid (2.0) expected 
for isotropic case. This is combined with a change of the 
shape of transverse section from a square of rounded cor¬ 
ners to a circle. Apparently, there appears to be some 
sort of a correlation between the tip radius and the num¬ 
ber of spiral arms: the probability of having larger num¬ 
ber of spiral arms increases with increasing tip radius. 
We note that the reason for observing v > 2 is proba¬ 
bly finite-size effect: the influence of the boundary of the 
simulation box on the dendrite shape exerted far from 
the tip region. 


D. Eutectic patterns 


Next, we investigate the formation of two-phase pat¬ 
terns on the surface of eutectic dendrites. As reported 
previously4® they include target-, single- and multi-arm 
spiraling motifs (Fig. 11). A closer inspection of the tip 
region reveals that several modes of pattern formation 
are possible. 

The target pattern appears via alternating nucleation 
of the two solid phases on top of each other [Fig. ©a)], 
confining the occurrence of this mode to larger undercool¬ 
ings. However, no nucleation is needed for the formation 
of single- and double spirals , in which the eutectic pat¬ 
tern originates from alternative occupation of the tip re- 
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FIG. 10. The effect of kinetic anisotropy on the shape of the 
two-phase dendrite, (a) Tip radius vs. kinetic anisotropy at 
7S£,o = 0.0295; (b) Exponent v vs. kinetic anisotropy. (v has 
been obtained by fitting the expression z = z ma x — \x\ v to the 
longitudinal cross section of the dendrite, where Zmax is the 
tip position, x is the distance from the axis of the dendrite.) 


gion by the two solid solution phases. In the latter cases 
this happens so that the a and /? phases remain spatially 
continuous: It appears that the alternating occupation 
of the tip region by the two phases may happen with¬ 
out nucleation. Having a growing a region at the tip, 
component B piles up ahead of the tip until fast growth 
initiated by the previously formed /3 phase leads to the 
formation of a /3 domain at the tip, which is thus con¬ 
nected to the rest of the (3 phase. This is followed by 
the the same process, but now for the a phase. Since in 
this way the formation of spiral arms (single or multiple) 
does not need nucleation, these mechanisms can prevail 
at small undercoolings. The growth modes with larger 
number of spirals (N aim > 3) become increasingly more 
complex, still displaying alternating phase appearance at 
the tip. Owing to the complexity of the process, it is 
difficult to decide whether nucleation plays a role here or 
not. The larger is the number of spiraling arms, lV arm , 
the more disordered is the tip region. As a result, larger 
numbers of point- and line defects travel down on the 
surface of the dendrite. This is especially so for the sim¬ 
ulations performed with flux noise for the concentration 
fields. 

The internal structure of eutectic dendrites with the 



FIG. 11. Eutectic patterns observed under nominally the 
same conditions (reference conditions as specified in Section 
III): (a) Target and, (b) single- and (c)-(f) multiarm spiral 
motifs. They differ in only the initialization of the random 
number generator used in creating the initial random eutectic 
pattern. 


target pattern, the single-, triple-, and five-arm spiral¬ 
ing modes are presented in Fig. |12| Apart from period¬ 
ically occurring defect-rich regions, the target pattern is 
composed of conical domains of the two phases, where 
the conical domains of the same phase are not connected 
spatially (i.e., alternating nucleation of the two phases 
is needed to create them). In contrast, the single- and 
multiple spiraling patterns observed on the surface are 
realized by single- or multiple helical structures forming 
in the volume. Owing to obvious geometrical constraints, 
the steepness of the spirals on the surface increases with 
the number of spiral arms. Remarkably, the longitudi¬ 
nal sections of the two-phase dendrites are fairly similar 
for all these modes, although weak systematic differences 
can be observed [Fig. [l2|]e)]. More characteristic are the 
front views and the transverse sections: The individual 
eutectic growth modes (number of spirals) can clearly be 
distinguished [see Fig. (Ha) -(d)]. It has been observed 
that in the appropriate parameter domain, the spiraling 
two-phase dendrites are quite robust. It is worth noting 
that in the experiments performed at low undercoolings 
the single-spiral mode has been exclusively observed 1 -. 
The multiple spiral arms seen in our simulations probably 
follow from the large relative undercooling we applied. 

Although the number of spiraling arms (-/V ar m) tends 
to increase with increasing solid-liquid interfacial free en¬ 
ergy (Fig. E- which in turn is reflected in increasing tip 
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Dendrite cross-sections for the target, single-, 

(a)-(c): one phase 
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FIG. 12. 

triple- and five-arm cases shown in Fig. 
x — y plane frontal view, the box cut at z = 3, 50 and 100, 
respectively, perpendicular to z axis. The eutectic struc¬ 
tures are connected inside the dendrite, (d) Helical structure 
formed by one of the phases. The cut surfaces were colored 
with blue for better contrast, (e) y — z plane cross-section. 


radius, the steady-state pattern depends also on the ini¬ 
tial random distribution of the two solid phases repre¬ 
senting the cumulative effect of the preceding composi¬ 
tional fluctuations: Different steady-state patterns were 
obtained starting from different (random) initial patterns 
(see Fig. 13). For example at 7 sl,o = 0.059, these pat¬ 
terns include the target pattern together with single- to 
fivefold spirals, which indicates the multiplicity of the 
possible steady-state solutions for the same physical con¬ 
ditions of which the random initialization chooses. 

This stochastic behavior is characterized by a peaked 
probability distribution (see central and bottom rows in 
Fig. 13), which shows increasing number of spiral arms 
with increasing solid-liquid interfacial free energy (and 
therefore tip radius). 

In order to check the validity of the assumption that 
such random initialization can reasonably represent the 
effect of compositional fluctuations, we have performed 
20 simulations starting from a single-phase crystal seed 
while adding flux noise to the EOMs for the concentra¬ 
tion fields as in Ref.^. These simulations differed in only 
the initialization of the random number generator. Such 
simulations are rather time consuming since growth from 
the seed to the steady-state dendrite has to be covered, 
and the adding of flux noise also slows down the simula¬ 
tion process. The corresponding probability distribution 
is compared to the distribution obtained from computa¬ 
tions with random initial two-phase patterns in Fig. |14| 
Considering the scattering of the results due the rela¬ 
tively small number of simulations limited by the avail¬ 
able computational power, a reasonable agreement can be 
seen between the two types of probability distributions. 
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FIG. 13. Pattern selection for two-phase dendrites: Top: 
Number of spiral arms vs. tip radius. Empty (green) cir¬ 
cles: results from random initial eutectic patterns at £4 = 0.3 
(180 simulations). Blue diamonds: same initialization, but 
e 4 varies between 0.3 and 0.05. Red squares: the most likely 
eutectic patterns from 20 random initializations (evaluated 
from the respective 20 empty circles). Bottom: Probability 
distributions (histograms) for steady-state eutectic patterns 
obtained from simulations started from 20 different random 
initial eutectic patterns. 


Apparently, in our case the thermal fluctuations decide 
which steady-state solutions are accessible for the system 
for a given set of operating parameters. These features 
closely resemble the behavior reported for Liesegang pat¬ 
terns, where the fluctuations determine, which of the 
competing modes (rings, single-, double-, tripl e- an d mul¬ 
tiple helices, or disordered patterns) is realized^- 13 . These 
findings raise the possibility that stochastic pattern se¬ 
lection is universal for certain multiarm spiral systems, 
in which diffusion plays a decisive role. Further work is 
needed, however, to clarify this issue. 


E. Is anisotropy needed for eutectic dendrites? 

One of the most interesting implications of the work 
of Akamatsu et a/.® is that unlike in the case of single¬ 
phase dendrites, where the presence of anisotropy is a 
precondition of dendritic growth®®®, it is expected that 
in the case of spiraling ternary eutectic dendrites den- 
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FIG. 14. Probability distributions (histograms) for the ac¬ 
cessible eutectic patterns differing in the number of spiral 
arms obtained on the basis of 20 simulations. Red: started 
from random initial eutectic patterns without noise. Yellow: 
started from a single-phase seed while applying flux noise in 
the EOMs for the concentration fields. The individual runs 
differed in the initialization of the random number generator. 
Note the similarity of the results obtained with and without 
flux noise for the concentration fields. 


dritic growth is possible without anisotropy. In order to 
clarify this issue further, we have performed simulations 
for £4 = 0. To minimize finite size effects that, in princi¬ 
ple, might contribute to the realization of the two-phase 
dendrites, we have made large cross-section simulations 
(on 512 x 512 x 200, 380 x 380 x 511, and 256 x 256 x 800 
grids) at the CCMS, IMR, Tohoku University. These sim¬ 
ulations were started from 25 to 100 randomly positioned 
single-phase seeds. Performing the simulation without 
anisotropy at v p = 0.1, the growing eutectic particles 
impinge upon each other, and after a transient period 
of frustrated eutectic patterns, the solid-liquid interface 
forms a 3D disordered cellular morphology with dynami¬ 
cally appearing and disappearing tips and ridges of rather 
flat end, covered by a disordered lamellar eutectic pattern 
with lamellae perpendicular to the solid-liquid interface 
(Fig. |15| ) ; patterns well known as eutectic colonies seen in 
ternary systems experiment all jEH and theoreticall y 36 * 77 ! 
The effect of pulling velocity and anisotropy are illus¬ 
trated in Fig. [lG| The growth forms in the isotropic sys¬ 
tems are shown in Figs. p~6|) a) and (b). In both cases 
dynamically changing cellular structures are observed, 
however, for the larger velocity the ’’lamellae” are close 
to perpendicular to the temperature gradient, indicat¬ 
ing a nearly ID eutectic growth expected at extreme 
high velocities^! a mechanism possibly emerging from 
the fact that at large undercoolings the applied model 
has a spinodal mechanism to form the solid phases upon 
each other. While in case (b) local spiraling may appear 
temporarily at the tips , well developed steady-state spi¬ 
raling dendrites cannot be seen. In these isotropic cases, 
the larger flat tips tend to undergo a tip-splitting/tip 
elimination phenomenon akin to the dynamic tip split¬ 


ting and cell elimination process shown by quantitative 
phase-field mode ling o f cellular growth in binary alloys 
of low anisotrop y 75 ^ 76 ! The respective simulations per¬ 
formed with kinetic anisotropy of £4 = 0.3 are displayed 
in Figs. 16 (c) and (d). While for the smaller velocity 
a dynamically changing cellular structure is observed, at 
the higher velocity [Fig. [l6](d)] , the peaks are far more 
regular than in the other cases and converge to a sta¬ 
ble dendritic configuration. Remarkably, we do not see 
side branches even in these large simulations. (We sus¬ 
pect that considerably larger simulation box is needed to 
capture this feature.) It appears that the spiraling eutec¬ 
tic pattern is considerably more definite with anisotropy 
than without anisotropy [cf. Figs. [I6|(a), (b), and (c) 
to Fig. |l6|(d)]. Note that these results, do not rule out 
the possibility that in the absence of anisotropy spiraling 
steady-state eutectic dendrites occur as excessively rare 
events. Large scale computations are needed to clarify 
this issue further. 


V. SUMMARY 

We have investigated the formation of ternary eutec¬ 
tic dendrites within the phase-field theory. It has been 
shown that a minimal ternary phase-field model of eu¬ 
tectic solidification suffices for reproducing the spiraling 
eutectic dendrites. We have mapped the domain of the 
Gibbs simplex, in which eutectic dendrites of ordered pat¬ 
terns (spiraling or target) appear. Increasing the pulling 
velocity, we see the following sequence of transitions be¬ 
tween morphologies/patterns: flat front lamellae —I eu¬ 
tectic colonies —I eutectic dendrites —I dendrites with 
target pattern —> partitionless dendrites (due to solute 
trapping for components A and B, but not for C) —> flat 
partitionless solid (due to full solute trapping). We have 
confirmed the assumptions of Akamatsu et alW^ that (i) 
the Jackson-Hunt scaling is followed by the spiraling eu¬ 
tectic pattern; and (ii) that the eutectic wavelength and 
the tip radius of the dendrite are proportional and of 
comparable magnitude. The geometrical shape of the 
eutectic dendrites and the tip radius behave analogously 
to their single-phase counterparts. Apparently, the un¬ 
derlying eutectic pattern has little influence on the shape 
of the dendrite. A number of quasi steady-state eutectic 
patterns have been observed under nominally the same 
conditions, including the target pattern, and single to 
multiple spirals, of which the fluctuations choose. This 
stochastic behavior resembles to that of helical Liesegang 
system^TS, raising the possibility of a universality class 
for systems displaying spiraling under diffusion control. 
It remains to be seen, however, whether this stochastic 
behavior is indeed general. In spiraling dendrites the a 
(and also the /3) phase domains are continuous through¬ 
out the solidified regions. It appears that the number of 
the spiral arms increases with increasing tip radius. The 
expectation that no anisotropy is needed for the forma¬ 
tion of two-phase dendrites 31 has not been born out by 
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FIG. 15. Formation of eutectic colonies in the model ternary eutectic system. Snapshots of a large scale simulation (380 x 
380 x 511 grid) performed at v p = 0.1 without anisotropy (taken at 10 4 ,6.2 x 10 s , 1.6 x 10 6 , and 4.0 x 10 6 time steps). Here 
solidification was started from 25 single-phase seeds. Note the competing cellular structures evolving after a transient period, 
and that the cellular structure dynamically changes with time. 


our simulations: Without anisotropy, the emerging cellu¬ 
lar structure does not display regular steady-state spiral¬ 
ing. It cannot be excluded, however, that in the isotropic 
systems spiraling eutectic dendrites are rare events that 
cannot be easily captured by simulations. 
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